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Abstract 

We use the linear a model to analyse the dynamics of a disoriented chiral 
condensate. For idealized boundary conditions appropriate to high energy 
collisions, the problem can be reduced to a one dimensional one. The evolu- 
tion of the chiral state is then that of a simple dynamical system and can be 
studied analytically. 
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1 INTRODUCTION 



The aim of this paper is to present a simple, analytical discussion of the time 
evolution of a disoriented chiral condensate (DCC) which could be produced 
in a high-energy collision 0-0. We shall extend our previous investigation 
H and make contact with recent works which have been triggered by the 
idea of a "quench", put forward in ref. ||. Specifically, we shall study the 
classical equations of motion of the linear a model, assuming random initial 
conditions for the fields and their derivatives. This randomness could result 
from a sudden cooling of a hot quark-gluon plasma, as proposed in [||], but 
we shall not insist on an interpretation based on initial thermal equilibrium. 
The evolution of the system will be described by the classical equations of 
motion of the a model. We neglect possible quantum effects (for a recent 
discussion of the decoherence of the radiation off colliding heavy ions see Q). 
The main motivation for using the linear, instead of the non-linear a model 
as in ||, is to allow for the possibility to start the evolution from a chirally 
symmetric random configuration, as in [|J. As we shall see, however, many 
features of our previous study survive. 

As in ref. || we adopt the idealization originaly proposed by Heisen- 
berg M: we assume that at time t — the whole energy of the collision is 
localized within an infinitesimally thin slab with infinite transverse extent. 
The symmetry of the problem then implies that the fields depend only on 
the invariant r = Vt 2 — x 2 , where x is the longitudinal coordinate. This 
is of course a dramatic simplification of the actual situation. In particular, 
our assumptions lead to unrealistic correlations both in the transverse plane 
and in rapidity. However, with this idealization the problem becomes 1+1 
dimensional, and can be treated analytically. This will allow us to identify 
clearly various regimes of the dynamics of chiral condensates in the presence 
of the longitudinal expansion which necessarily occurs in case of high en- 
ergy collisions. We believe that such analytical studies are complementary 
to numerical simulations [^|, |7|, aiming at getting a more realistic picture, 
and they may prove useful in understanding the results obtained in such 
simulations. 

Although we use Heisenberg's idealization to justify working in 1+1 di- 
mensions, we restrict our discussion to proper times r > r . An extension 
of the model down to r = leads to unphysical results and is completely 
unjustified. Heisenberg || has pointed out that in a wide class of non-linear 
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theories with gradient interactions, the fields have to vanish like r in the 
limit r — > 0. It seems more likely, however, that the fields (7?, a) do not tend 
to any well defined limit as r — > but behave chaotically (in the loose sense 
of the word). In any case, as r — > , the relevant degrees of freedom are 
presumably no longer hadronic and, besides, the effective theory is certainly 
more complicated than the linear a model. 

Thus, we shall assume that starting from some finite time r , we can 
describe the state of matter produced in a collision by the linear a model. 
Our results will turn out to be little dependent of the specific value of this 
time To, as long as it is chosen so that at tq the system is in the chirally 
symmetric configuration. The values of the fields at this time are taken as 
random variables. In fact, as we shall show in section 2, the equations of 
motion are, for short times, equivalent to the conservation of the isospin 
vector and axial vector currents. The bulk of the dynamics will be entirely 
determined by the values of these currents at t . In section 3, we shall identify 
several regimes of the evolution of our chiral condensate. We shall see that 
the problem reduces to that of finding the motion of a fictitious particle in 
a time dependent potential and subject to a time dependent friction. The 
paper ends in section 4 by a discussion. 

2 FORMULATION OF THE MODEL 

Our system is described by the Lagrangian of the linear a model, with the 
standard chiral symmetry breaking term: 

L = l -[{daf + («9tt) 2 ] - -y + Tf 2 - I) 2 + Ha. (1) 

The fields a and 7? are dimensionless and the parameters A and H have 
the dimension of [mass] 2 . We have scaled out a dimensionless factor S 1 / 2 , 
which would normally occur in the expression of the Lagrangian for the 3- 
dimensional model M. The parameters of the model may be related to the 
masses of the a and pion field, respectively: m 2 ~ 2A and m 2 = H. For the 
sake of the argument we shall proceed as if one had m a 3> m n . We shall 
return to this point in the discussion. 
The equations of motion read 
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T -\raJ=-X(cx 2 + n 2 -l)a + H, (2) 

and 



T T7T 



A(cr + 7T — 1)7?, (3) 



where the prime denotes differentiation with respect to r. From (|3p one easily 

gets 

7T X 7?' = (4) 

r 

where a is an integration constant. This is a direct consequence of the conser- 
vation of the isovector current together with the assumption that the fields 
depend on r only. Hence the component of 7? along a vanishes: 

7T a = 0. (5) 

A further consequence of the equations of motion is the equation 

T~ l [r(7ra' - an ')]' = Hit, (6) 

which is just the statement of partial conservation of the isoaxial- vector cur- 
rent (PCAC). 

For small enough r, the effect of the symmetry breaking term Htt can be 
neglected. The isoaxial-vector current is then conserved, and one can write, 
in analogy to (f|) 

7?Cr' — (J7T ' = — , (7) 

r 

where b is another integration constant. When eqs. (§) and (Q) are satisfied, 
the component of 7? along c = a x b equals a, up to the constant factor a/b. 
The forms of the conservation laws (^) and (0) reveal important features of 
our model. Because of the symmetry of the fields, and the classical dynamics, 
in the limit of vanishing pion mass the isospin vector and axial vector currents 
keep fixed orientations in isospace as the system evolves. Of course, any 
orientation is equally probable. But once fixed in the initial conditions, it 
remains constant. 
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The equations (||) and (0) are identical to those obtained in the non- 
linear a model 0, and they are independent of A. The non-linear a model 
is obtained in the limit A — > oo. In this limit, the motion takes place on 
the hypersphere a 2 + if 2 = 1 and is entirely determined by the conservation 
laws for the isospin currents. When A is large, but not infinite, one can still 
separate the general motion into a motion on the sphere, which one may 
loosely associate with the pion degrees of freedom, and a motion normal to 
the sphere which involves the a degree of freedom. As long as the parameter 
H can be neglected, there is no scale controlling the motion on the sphere, 
while A -1 / 2 controls the time dependence of the oscillations normal to the 
sphere. When H can no longer be ignored, it controls the time behaviour of 
pionic excitations. 

In order to analyse the evolution of the chiral state, we shall use the 
following parametrisation which takes into account the symmetries of the 
problem: 

Tib = r sin 9 

7r c = r cos 9 sin uj 

a = rcos9cosu. (8) 

Note that a 2 + tt 2 = r 2 . Our problem has now been reduced to a simple 
dynamical problem for a point particle with one radial degree of freedom, r, 
correesponding to the a excitations, and two angular degrees of freedom 9 
and u, corresponding to the pionic excitations. As will become clear later, 
the evolution of the system depends upon the initial conditions mostly via 
the constants a and b (while the detailed initial values of r, 9, uj and their 
derivatives are not essential). The two constants a and h determine the 
strength of the isospin vector and axial vector currents (note that due to 
a different choice of variables, the present a and b are twice as large as 
those introduced in ref.0). They are random variables whose probability 
distribution will be calculated in the last section. 

As a final remark in this section, we note that when eqs. (|j) and (|7|) are 
satisfied, the angle uj is constant and given by 

a 

uj = arctan — . (9) 
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Thus the motion is planar in the 3-space (a, lib, tt c ). This plane becomes 
almost identical to (TTb, er) when b 3> a. We shall find convenient, at some 
point in the forthcoming discussion, to work in this particular limit. 

3 DCC AS A SIMPLE DYNAMICAL SYS- 
TEM 

The initial time r is supposed to be such that the mass term oc H can be 
ignored at the beginning of the evolution. Then the initial motion is governed 
by the two conservation laws for the isospin currents, i.e. eqs. (^p and (0). 
These translate into the following equation for the "angular momentum" of 
the fictitious point particle associated with our system: 

r 2 9' = -, (10) 
r 

The equation for the radial motion is obtained by neglecting H in the equa- 
tions of motion (^)-(|3|), multiplying the equation for a (ti^) by cos 6 (sin#), 
and adding the resulting equations. We get 

r" + lrf = jf-- Ar(r 2 - 1). (11) 

We have used eq. (p~0|) to eliminate 8', and also the fact that uj is constant 
when the isospin currents are strictly conserved. The second term on the 
left-hand side may be interpreted as a time dependent friction. Thus, the 
fictitious particle experiences damped oscillations. The amplitude of these 
oscillations can be initiallly quite important if r(r ) is small. In that case, 
the repulsive force represented by the first term on the right-hand side of 
( |TTD provides to the particle a large velocity. However, the oscillations are 
rapidly damped, while the repulsive force itself decays fast with time. 

Further insight on the radial motion can be gained by considering the 
time variation of the mechanical energy of the fictitious particle. This is 
given by 



This derivative is particularly large, in absolute value, during the very early 
stage of the evolution due to the presence of the term ~ t~ 3 . The whole 
damping is of kinematical origin; it reflects the decrease of the field energy 
density as the system expands. One observes that the influence of the angular 
motion on the radial motion quickly disappears. The energy in the angular 
motion is controlled by k. It becomes of the same order of magnitude as the 
potential energy oc A when r ~ «/\/2A, and is negligible at later times. 
For not too small r it is certainly a good approximation to linearize 



the expression appearing on the right-hand side of (|TT[) . There results an 
inhomogeneous Bessel equation for r — 1. The solution takes a particularly 
suggestive form for v2Ar ^> max[ft 2 , 1]: 

C 

r«H r cos(v / 2Ar + 77), (13) 



where C and 77 are integration constants fixed by the initial conditions. Thus, 
because of the damping caused by the expansion, the system is driven towards 
the orbit r = 1 [strictly speaking, for H 7^ 0, the orbit would rather be at 
r = 1 + 0(H/X)]. 

When r ~ 1, the integration of eq. (|10|) yields 9 « «ln(r/ro) and one 
recovers the solution of the non-linear a model presented in ref. ||. In fact 
the approach to the orbit is fairly slow, and a change of regime will in most 
cases take place before the orbit is actually reached. The oscillations in the 
radial motion have a large frequency, oc v2A. Since the angular motion, 
controlled by k, is generically slower than the radial one, one can, with a 
good accuracy, replace r 2 in eq. (|10D by its expectation value over a few 



periods of the radial motion. This expectation value reaches a value 1 in 
a short time. Therefore, the result 9 m k\ti{t/tq) can be taken as a good 
approximation from nearly the beginning of the evolution. 

A change of regime is expected to occur when the mass term oc H can 
no longer be ignored. Let us, therefore, discuss in more detail the quality of 
the approximation H = 0. Integrating eq. (|6]) for TTf, one has 

r (717,0"' — an' b ) = b + H J rr sin 9 dr. (14) 

In estimating the integral above it is important to notice that the integrand 
is oscillating, and it does so quite rapidly when k is large. In order to 
find a rough estimate, we replace again r by 1, its average value, so that 
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9 = k In (t/tq) (we assume 0(tq) = to avoid the proliferation of inessential 



integration constants). The integral in (14) is then easily evaluated. For large 
k, it equals (r 2 //t) cos (ft In (t/t ). Therefore, the neglect of H is justified as 
long as 



T 



bn 



« T (15) 

It remains to determine the conditions ensuring that u = const., because the 
planarity of the motion in the 3-space (a, 7tb, 7i c ) is another characteristic fea- 
ture of the early, chirally symmetric regime. Using the polar parametrisation 
(|[) and the equations (f|) and (|]) one finds after some algebra 



cos 9 sin (lu — uq) = — [cos 9 sin uj I rr sin 9dr — sin 9 rr cos 9 sin ujdr] 

(16) 

where ujq = arctan (a/b), the value uj takes when H = 0. On the right-hand 
side set uj = ojq + 5u and expand neglecting terms of 0[(5o;) 2 ]. Furthermore, 
set again r = 1 and 9 = /tin (r/r ) to get 

Ht 2 bH r 
cos 9 6 <oj ~ ^-[a — 6cos 2 ^] r-sin^ / r cos 95ojdr (17) 

One can see that, as long as cos# is not small, \8u\ is of order Hr 2 /n 2 . 
However, when cos# = 0, the left-hand side of flTTD vanishes identically, and 
a perturbative (in H) calculation of 8uj is impossible. However, one easily sees 
that the right-hand side of (O) vanishes when the time average (5uj) rs —a/b. 
We conclude that the angle uj oscillates with an amplitude of order a/b, i.e. 
independent of H. We shall assume herefrom that 

6>a (18) 

The constraint fllTf ) makes the problem technically simpler, since uj remains 
small and the motion on the orbit is quasi-planar at any time. We do not 
think that limiting the discussion to this case is a severe restriction. The 
general case is just more cumbersome, without being really more instructive. 
With (|T8|) we gain considerably in clarity. Furthermore, only for large enough 
b the problem is really interesting, a point that will become more and more 
clear as the discussion will develop. Finally, as will be shown in the last 
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section, it is improbable that both b and a are large. When (|I8| ) holds, eq. 
(|15D can be rewritten as follows 



r « (19) 



The presence of H in the equations of motion is an inessential complica- 
tion as long as ( |19|) is true. We shall see that a dramatic change of regime 
occurs when r « b/y/H. 

Consider first what happens at large times. To this end we use eqs. @) 
and Working in the polar parametrization (P), with r = 1, one easily 
convinces oneself that 9, uo — > like some inverse power of r when r — > oo, 
modulo possible oscillations. For small enough 9, u eqs. (0) read 



t-^tO')' = -HO 
t -\tu')' = -Hu, 

while eq. (|]) becomes a constraint on the Wronskian W(u,9): 

W(u,9) = -. 



(20) 



(21) 



Recall that the operator appearing on the left-hand side of eqs. ( BD[ ) is the 
form taken by d 2 acting on a field depending on r alone. Furthermore, in 
the regime that we are now considering, 9 ~ 717, and u ~ 7r c . Hence, eqs. 
([20"D describe the free propagation of massive pions, the mass being if 2, as 
expected. Eqs. (|20D become Bessel equations when the variable V^ffr is 
used. One has 



ffr) 



T cos (VHt + C) 



D(VHt) 



(22) 



for ^> 1, where -D and C are integration constants. Notice, that the 

amplitudes and phases of 9 and u are related by the constraint (|2l|). Eqs. 
( |22|) hold asymptotically, independently of the values taken by the constants 
a and b. 
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In order to determine the constant D one needs to match the two ex- 
treme regimes discussed so far, which is not an easy task if attempted in full 
generality but which simplifies considerably when a ^b. 

By using the approximation already discussed in which one replaces r by 
its expectation value, and rr' by zero, and assuming that uo remains small 
enough, we can transform eq. (^|) for -n^ into 

0" + + # sin = 0. (23) 

T 

This is the equation of a pendulum subject to a time dependent friction. The 



asymptotic solution, for r — > oo, has already been written in eq. ( 22|) and 
describes an oscillatory motion. The initial value of is not very important. 
However, for large b, the derivative 0' = b/r remains large, generically, at 
the time when the expectation value of r becomes close to unity. As long 
as the angular velocity is large enough, the motion is circular. Our purpose, 
now, is to discuss the transition from the circular to oscillatory motion of 
this simple dynamical system. In particular, we shall determine the time 
when the transition takes place and the amplitude of final oscillations (the 
parameter D in eqs. (|2~2"D). 



Eq. fl23|) implies that 

Tr E{T) = Tr [ \ {6 ' )2 + H(1 " COS9)] = -\ {e ' f < °- (24) 
Thus the "energy" E{r) of the pendulum is a strictly decreasing function of 
time. As long as the total energy is much larger than the potential energy 
H(l — cos0) , the motion is circular. The transition occurs, when E{r) "hits" 
the potential. Let us solve eq. (0), treating this potential as a perturbation. 
One finds 



0^61n(r/r o ), (25) 

and 

0' » - + F^cos[6 In (t/t )]. (26) 
r o 

Since the slope of E(r) decreases rapidly with r, for large b the potential is 
"hit" near its top at time r satisfying cos [bin (t/to)] ~ —1. The transition 
time is estimated by setting 0' = and is found to be 
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Ttrans « ~f=, b > 1. (27) 
V 

As one might have expected, the transition occurs when the naively estimated 
kinetic energy becomes comparable to the potential energy. The amplitude 
of oscillation that results is, of course bounded by it. Actually, it is close to 
7r just after the transition, precisely because the potential is "hit" near its 
top. Linearizing the potential one gets the solution fl22[) . Keeping in mind 
the remark just made about the amplitude, and using fl27D one finds 

const. Vb, 6>1, (28) 

where the constant is of order unity. The amplitude of the pion field is oc y/b 
and, consequently, the height of the rapidity plateau is oc b. This is similar 
to the result found in |J in the framework of the non-linear a model. As was 
the case there, the parameter b (« k when a <C b) is a measure of the energy 
released in the decay of the condensate. 



4 DISCUSSION 

We have identified several stages in the evolution of the chiral condensate. 
First there is a short phase in which angular motion and radial motion are 
strongly coupled. This phase lasts for a time ~ b/\/2X. Next, the ficti- 
tious particle rotates slowly, while oscillating rapidly about the equilibrium 
radial position r « 1. Because of the expansion of the physical system, the 
energy of the fictitious particle is damped. At some point there is a transi- 
tion from circular to oscillatory motion, which actually corresponds to free 
propagation of final state pions. The transition takes place when the mass 
of pion excitation can no longer be neglected, and it occurs typically when 
t ~ b/VH. 

The pion field has a random orientation in isospace and therefore the 
fraction / of neutral pions is distributed according to the law 

dP(f) = ^j=df (29) 

written first explicitly in ref. || (but obtained earlier, and independently, by 
Bjorken ||). At all times the system is subject to "friction". The latter is of 
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purely kinematical origin and reflects the boost invariance imposed on our 
solution. There is no true energy dissipation. Simply the expansion causes 
the field energy in a covolume to decrease. 

The two time scales T\ = (2A)~5 and r 2 = (H)~2 y if estimated using the 
phenomenological values of the pion and a masses, differ only by a factor 
of 4. This is presumably not sufficient to guarantee a clean separation of 
regimes, especially when b is not large enough (and may be related to dif- 
ficulties in producing large domains of misaligned vacuum reported in ref. 
[0]). Although it seems that the picture becomes sharper when the conden- 
sate carries more energy, there is a price to be paid for that: as will be shown 
in a moment, the probability that b is large falls exponentially. 

Let the fields (7?, a) and their (proper) time derivatives be Gaussian ran- 
dom variables with variances o\ and a^, respectively, at some initial time To 
belonging to the chiral symmetric regime (cf. ref. H). With this assumption 
let us calculate the probability P(n, uj) to find a given value of k = \/a 2 + b 2 
and u = arctan(a/fe). Using the polar parametrization of the fields in the 
randomly oriented hyperplane _L a, together with the constraint k = T§r 2 B' 
(cf. (0)), one finds, after a straightforward calculation: 

dKaio 2{7r<Jia2T ) z VeriO^To 

where K\(z) denotes the second modified Bessel function, K\{z) ~ <Jir/2ze~ z 
for z — > 00. The probability distribution (|3~0"D falls exponentially at large k. 
This result is unlikely to be of much phenomenological relevance, since it 
is obtained in a 1 + 1 dimensional model. The reduction to 1 + 1 dimen- 
sions means full translational symmetry in transverse coordinates and, in 
particular, coherence over large transverse distance. The probability that 
such coherence does occur cannot be estimated within the framework of the 
model. Eq. ([30]) represents nevertheless a significant conceptual progress 
with respect to the discussion of refs. -[||: the present approach, which 
has been partly inspired by ref. [f|], enables one to determine the absolute 
probability of observing a disoriented chiral condensate with given global 
characteristics. 

The initial randomness of the fields (7?, a) , and their derivatives, does not 
necessarily reflect thermalization occuring at the early stage of the collision 
process. Actually, the most interesting aspect of the quench model put for- 
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ward, in this context, by Rajagopal and Wilczek is the out of equilibrium 
evolution of the system. In discussing this evolution one does not need to 
talk about temperature at all. Since the seminal work of Landau [O, the 
study of complicated nuclear collisions uses the concept of (local) thermal 
equilibrium as the central piece of the paradigm. If we were to assume that 
local thermal equilibrium is maintained throughout the time development, 
we would be led to use hydrodynamical evolution equations. The resulting 
evolution would be very different from that obtained here. It is unclear, at 
this moment, whether the disoriented chiral condensate can be really pro- 
duced in high energy collisions. However, the theoretical investigation of this 
as yet undiscovered phenomenon has already turned out to be useful, since 
it has been a good pretext for trying to figure out non-equilibrium aspects 
of high-energy nuclear collisions. 
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